Stretching of vortex lines and generation of vorticity in the three-dimensional 

complex Ginzburg-Landau Equation 
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The dynamics of curved vortex filaments is studied analytically and numerically in the framework of 
a three-dimensional complex Ginzburg-Landau equation (CGLE). It is proved that a straight vortex 
line is unstable with respect to spontaneous stretching and bending in a certain range of parameters 
of the CGLE, resulting in formation of persistent entangled vortex configurations. The analysis 
shows that the standard approach relating the velocity of the filament with the local curvature, is 
insufficient to describe the instability and self-generation of vorticity. 
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The complex Ginzburg-Landau equation (CGLE), de- 
rived some 20 years ago by Newell and Kuramoto ||] 
has become a paradigm model for a qualitative descrip- 
tion of weakly- nonlinear oscillatory media (see for review 
||). Under appropriate scaling of the physical variables, 
the equation assumes the universal form 



d t A = A - (1 + ic)\A\ 2 A + (1 + ib)AA, 



(1) 



where A is the complex amplitude, b and c are real pa- 
rameters, and A = <9 2 + <9 2 + <9 2 is a three-dimensional 
Laplace operator. Although the equation is formally 
valid only at the threshold of a supercritical Hopf bi- 
furcation, it has been found that the CGLE often repro- 
duces qualitatively correct phenomenology over a much 
wider range of the parameters. As a results, the predic- 
tions drawn from the analysis of the CGLE (mostly in 
one and two spatial dimensions, see e.g. ^-§]) were re- 
cently successfully confirmed by experiments in optical 
and chemical systems [0J|]. Moreover, some results ob- 
tained from the CGLE (for example, symmetry breaking 
of spiral pairs) was instructive for interpretation of ex- 
periments in far more complicated systems of chemical 
waves H and colonies of amoebae jlC],|llJ . 

Recently, the dynamics of three-dimensional (3D) 
vortex filaments has attracted substantial attention 
||,0-[l7). In the context of the 3D CGLE, Gabbay 
Ott and Guzdar |l4|] applied a generalization of Keener's 
method for a scroll vortex in reaction-diffusion systems 
p8[ . They derived that a vortex ring of radius R col- 
lapses in finite time over the entire range of parameters 
(6, c) of Eq. (0) according to the following evolution law 



dR 
~dt 
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This result generalizes Keener's ansatz by including the 
curvature- induced shift of the filament's wavenumber and 
corrects the erroneous estimate of Ref . [O . Thereby, as 



follows from Eq. (g), vorticity initially existing in the 
system will always decay (if, of course, there is no bulk 
instability of the waves emitted by the vortex filament), 
and under no condition can the vortex filament expand. 
The persistence of vortex filaments, known under some 
conditions for reaction-diffusion systems, was attributed 
in the context of the CGLE to the effect of higher-order 
corrections to the CGLE. 

In this Letter we show that under very general con- 
ditions and in an extensive part of the parameter space 
vortex filaments expand and result in persistent vortex 
configurations even if there is no bulk instability of emit- 
ted waves. The condition for the expansion of the vortex 
filament is simply b > b c (c) 3> 1. In this limit Eq. (Q) 
is not valid, because formally higher-order corrections, 
omitted in Eq. (||), cause severe instability of the fil- 
ament and persistent stretching. This instability is a 
three-dimensional manifestation of the two-dimensional 
core instability of spiral waves (called in Ref. |^| accel- 
eration instability). It originates from breakdown of the 
Galilean invariance of the CGLE for any e = 1/i ^ 0, 
causing spontaneous acceleration of the spiral waves for 
e < e c [||- Although in 2D the instability is relatively 
weak (the growth rate is proportional to e) , the situation 
is different in three dimensions. As we will show, the 
bending of the filament greatly enhances the instability, 
resulting in formation, after some transient, of a highly 
entangled and dense vortex configuration. The condition 
b 3> 1 is readily fullfiled for many physical and chemical 
systems. For example, in the context of nonlinear optics, 
where the CGLE can be derived from the Maxwell-Bloch 



equation in the good cavity limit 13 1, this parameter is 
very small: e = l/b ~ 10~ 4 — 10^^ For an oscillating 
chemical reaction the diffusion rates of various compo- 
nents can be varied over a wide range by adding extra 
chemicals. 

As a test for instability, we consider the dynamics of a 
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weakly curved vortex filament. We apply the generaliza- 
tion of the method of Ref. || for the case of 3D vortices, 
and make perturbations near the 2D spiral wave solution 
to the CGLE. For convenience, we redefine r — > r/Vb. 
Then Eq. (fy) assumes the form 



d t A = A - (1 + ic)\A\ 2 A + (e + i)AA. 



(3) 



In the following discussion, we assume < e <C 1 to be 
a small parameter. Our objective is to relate the accel- 
eration of the vortex filament <9(V with the velocity v 
and local curvature k of the filament. We obtain that for 
e< 1 the equation of filament motion can be written as 



d t v + K[ev + (I + £ 2 )kN] = 0, 



(4) 



where N is the unit vector pointing toward the center 
of curvature, and B is the unit vector perpendicular 
to the filament and N. Velocity and acceleration have 
correspondingly two components, v = (vn, Vb), <9jV = 
(dtV]\[,dtVB)- K is the friction tensor satisfying Kn = 
K22,Ki2 = — K21. Dropping the first term in Eq.([|), 
we reproduce the result of Ref. jlij for the collapse rate 
v jy = —(1 + e 2 )n/e (because K = l/R)- However, if 
Kn < 0, which we will show is the case for e < e c (c), 
the acceleration has a significant effect and leads to an 
ultimate expansion of the vortex filaments and the per- 
sistence of entangled vortex configurations. In principle, 
the numerical value of e c (c) does not need to be very 
small. For example, for c = 0.5 we obtained e c ~ 0.165, 
which corresponds to b rs 6. 

In order to develop the perturbation theory for a 
weakly-curved vortex filament in 3D, we begin with the 
stationary one-armed isolated spiral solution to Eqs. 
(plph, which is of the form 



A (r, 9) = F(r) expi[ut ±9 + xjj(r)} , 



(5) 



where (r, 9) are polar coordinates. The real functions F 
and ip have the following asymptotic behavior F(r) — > 
y/l — (k 2 , ip'(r) — ► k for r — > 00 and F(r) ~ r , ip'(r) ~ r 
for r — * 0. The wavenumber k of the waves emitted by 
the spiral is determined uniquely for given e,c |fj"9fl . For 
e = one has a type of Galilean invariance and then, in 
addition to the stationary spiral, there exists a family of 
spirals moving with arbitrary constant velocity v jq] , 



A(r, t) = F(r') exp i[uj't + 9 + tp(r') - 



(6) 



where r' = r + vi, u>' = u> + v 2 /4, and the functions 
F,ip are those of Eq. (||). (This invariance holds for any 
stationary solution) . For e ^ the diffusion term ~ e AA 
destroys the family and leads to acceleration/deceleration 
of the spiral proportional to ev, depending on the value 
of e. For e < e c (c) the spiral is unstable with respect to 
spontaneous acceleration since Kn <0|. 



Let us now consider the dynamics of an almost straight 
vortex line in the filament-based coordinate system (see 
for details The position in space X is represented 

by local coordinates s, x, y, where s is the arclength along 
the filament, and X = R(s) -KtN(s) + j/B(s), where R is 
the coordinate of the filament. In this basis the weakly 
curved filament moving with velocity v can be written in 
the form 



A(r,t) = F(r')exp 
S kg* x 



i (u't + 9 + ip(r') ~ 
+ W{r',6,s). 



(7) 



Here r' = r + vt, W is the perturbation to the spiral 
solution which we require to be small, and 5k x is the 
correction to the wavenumber due to curvature of the 
filament. In principle, this correction can be absorbed 
in W but it is convenient to retain this form since we 
can cancel part of the perturbation exactly by adjust- 
ing 5k x jl4j . The perturbation procedure to derive Eq. 
([|) is practically identical to that of Refs. [Q^]. Sub- 
stituting the ansatz Eq. (0) into the CGLE, and as- 
suming -C 1, one obtains from Eq. ([}|) in first 
order in e an inhomogeneous linear equation for the 
correction W, which is of the form LW = H, where 
LW = {l+iu-2(l+ic)\A \ 2 )W-(l+ic)A 2 l W*+iA ± W, 
A± = d 2 + d~. In first order the W-independent part H 
has two terms, H = H tr + H rot , where H tr contributes 
to the displacement of the vortex filament, and H rot con- 
tributes to the shift of the vortex frequency. We are 
interesting only in H tr (the contribution H rot is not re- 
sponsible for the instability), which is of the form: 



-r'dtv— + ievVAo + (e + - i8k x )S7A . 



(8) 



The derived system of equations is very close to that con- 
sidered in the Ref. || in the context of the acceleration 
instability of a spiral wave in 2D. The acceleration can 
be obtained as a result of the solvability condition requir- 
ing that W be regular at the core of the spiral and does 
not diverge exponentially at large r. (Slower, power-law, 
divergence of W is permitted). To impose the solvabil- 
ity condition a specific numerical procedure is required 
(see for details |§,||). Eq. (||) contains an extra term 
(e + i)(n — iSk x )d x Ao with respect to that considered in 
H, which originates from the Laplace operator in a lo- 
cal basis (in the limit of small curvature k and torsion 
t): A = — nd x + d 2 + d~ + d 2 + .... However, additional 
numerical integration is not necessary to account for the 
effect of this term, because, as was found in Ref. JM) , 
this term can be canceled exactly to first order in k by 
a proper choice of 5k x = —en. After this modification, 
H tr assumes the form H tr = — r'Stv^ 2 - +icvVAq, where 
vn = vn + (1 + e 2 )re/e and vb = vb- It is easy to see that 
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it now coincides with the corresponding function consid- 
ered in Ref. ||. Following the lines of analysis of Ref. 
||, we observe that the relation between dtv and v is 
identical to that considered in Q , and is of the form 

9 t v + eKv = 0. (9) 

Thus we obtain Eq. (^). The coefficients of the friction 
tensor coincide with those calculated in Ref. Q us- 
ing numerical matching, and for e < e c it was found that 
Kn < 0, which would guarantee the instability with re- 
spect to spontaneous acceleration of the vortex filament. 
Moreover, since in 3D the direction of motion of the fila- 
ment in general varies along the filament, this instability 
results in stretching of the vortex. 

In the case of the acceleration instability, K\\ < and 
the curvature has a strong destabilizing effect. Indeed, 
for an almost straight vortex parallel to the axis z, we can 
parameterize the position along the filament by the z co- 
ordiante: (Xo(z),Yq(z)). Since in this limit the arclength 
s is close to z, the curvature correction to the velocity kN 
is simply kN = — (<9 2 A , d 2 Yb). After simple algebra one 
derives the following relation for the growth rate X(k) 
of linear perturbation Xq{z),Yq{z) ~ exp[A(&)t + ikz]: 
A 2 + (Ku + iK 12 )(e\ + k 2 ) = 0. For k > e we ob- 
tain A ps y/—(Kn + iK\2)k 3> e. Thus, for finite k the 
growth rate X(k) may significantly exceed the increment 
of the acceleration instability in 2D (corresponding to 
k = 0): A = —e(Kix + %K\i). We can expect that, as 
a result of such an instability, highly-curved vortex fila- 
ments will be formed. Therefore, the above considered 
" small-curvature" approximation can be valid only for fi- 
nite time. Moreover, it is naive to expect a saturation of 
this instability in a steady-state configuration with finite 
curvature. In contrast, we suggest that frequently re- 
connection of various parts of the filaments, formation of 
vortex rings etc, will result in persistent spatio-temporal 
dynamics of a highly- entangled vortex state. 

In order to follow further development of the instabil- 
ity we performed numerical simulations of the 3D CGLE. 
We studied a system of 50 3 dimensionless units of Eq. (||) 
with no-flux boundary conditions. The numerical solu- 
tion was performed on an R10000 SGI workstation by an 
implicit Crank-Nicholson algorithm. The number of grid 
points was 100 3 . We performed simulations for e = 0.02 
in the parameter regime away from amplitude turbulence 
in 2D H . As an initial condition we used a straight vor- 
tex line perturbed by small noise. The results of the 
simulations are shown in Fig. 1. We see that the length 
of the vortex line grows. The dynamics seems to be very 
rapidly varying in time, and the line intersects itself many 
times forming numerous vortex rings. The long-time dy- 
namics show, however, some kind of saturation when a 
highly-entangled vortex state is formed and the length of 
the line cannot grow anymore due to repulsive interaction 
between closely packed line segments. The dependence 



of the line length on time is shown in Fig. 2. As a mea- 
sure of the filament length L we used the following quan- 
tity: L w Sq 1 J <d(A — \A(x,y, z)\)dxdydz, where 6 is 
a step function: Q(x) = 1 if x > and 9 = otherwise. 
Aq = 0.1 was used as a threshold value to identify the 
vortex. Sq is a constant determined from the condition 
that for the straight line the above integral coincides with 
the actual length. We can identify two distinct stages of 
the dynamics: first, fast growth of the length; second, 
oscillations of the line's length around some mean value. 
Remarkably, we observed an increase in the amplitude 
of the oscillation moving though negative c values in the 
regime of spatio-temporal intermittency in the 2D CGLE. 
It is plausible that the symmetry breaking mechanism re- 
sponsible for persistent intermittent behavior in the 2D 
system still has some importance in 3D. 

Persistent entangled vortex configurations are known 
from numerical simulations of excitable reaction-diffusion 
systems |ll|[l(|. Recently, elongation of a vortex fila- 
ment was also related with the meandering instability of 
a vortex core in 2D f20|| . However, there is also a signif- 
icant difference between these two phenomena. For the 
reaction-diffusion systems the meandering instability is 
typically supercritical, and may be saturated at a rela- 
tively small amplitude of oscillations. In application to 
3D systems, one expects that highly-curved vortex fila- 
ments will not develop, and for this reason the intersec- 
tions are very seldom. This picture is consistent with 
the numerical simulations. Consequently, one may ex- 
pect that the analysis performed on the basis of a small- 
curvature approximation will be valid for a very long time 
and can capture the actual dynamics of the filament. In 
the CGLE the 2D instability is subcritical and the veloc- 
ity of the spiral in an infinite system can grow arbitrarily 
until additional topological defects nucleate in the wake 
of the accelerating spiral. In the 3D situation we have 
shown that curvature of the filament is even a destabi- 
lizing factor. Our numerical simulations show that the 
filament does not approach any steady-state. In contrast, 
it shows "violent" intermittent behavior, with numerous 
reconnections and splitting of vortex rings. As a result, 
the approximate equation of motion (^) must be consid- 
ered only as a test for instability rather than long-time 
evolution. 

In conclusion, we have derived an equation of motion 
for the vortex filament in the CGLE. We have found that 
in a wide range of parameters of the CGLE the vortex 
filament is unstable with respect to spontaneous stretch- 
ing, resulting in the formation of persistent entangled 
vortex configurations. This emphasizes the deficiency of 
previous approaches relating local filament velocity to 
local curvature. Our result could be verified in exper- 
iments with autocatalytic chemical reactions in gels in 
the regime of oscillatory instability. Also, the limit of a 
large cross-diffusion ratio can probably be achieved by 
doping with additional chemicals, thus changing the rel- 



3 



ative mobility of one chemical species with respect to 
another. Recently, the amplitude equation governing the 
dynamics of an elastic rod was derived |2l| . We note that 
our instability can be formally interpreted as a dynamics 
of a thin rod with negative elasticity. We also speculate 
that our results are relevant for inviscid hydrodyanmics. 
In the limit of b, c — > oo, Eq. (|l|) reduces to the defo- 
cusing nonlinear Schrodinger equation (NSE), which is a 
paradigm model for compressible inviscid hydrodynam- 
ics. Although the vortex lines are stable in the framework 
of the NSE, the corrections arising from the CGLE cause 
the destabilization and self-generation of vorticity. 
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FIG. 1. Instability of a straight vortex filament. 3D iso- 
surfaces of \A{x,y,z)\ — 0.1 for e = 0.02, c = —0.03, shown 
at four times: 50 (a), 150 (b), 250 (c), 500(d). 
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FIG. 2. The dependence of filament length L on time. 
Solid line corresponds to £ = 0.02, c = —0.03; dashed line 
corresponds to e = 0.02, c = —0.5. 
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